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Abstract. In ultracold atoms settings, inelastic light scattering is a preeminent 
technique to reveal static and dynamic properties at nonzero momentum. In this work, 
we investigate an array of one-dimensional trapped Bose gases, by measuring both the 
energy and the momentum imparted to the system via light scattering experiments. 
The measurements are performed in the weak perturbation regime, where these two 
quantities - the energy and momentum transferred - are expected to be related to 
the dynamical structure factor of the system. We discuss this relation, with special 
attention to the role of in-trap dynamics on the transferred momentum. 


Submitted to; New J. Phys. 



CONTENTS 


2 


Contents 

1 Introduction [2] 

2 Energy and momentum transferred to an array of ID gases [3] 

2.1 Experimental setup. |3] 

2.2 Results. H] 

2.3 Comparison with the response of a 3D condensate. |6] 

3 Effect of in-trap dynamics on the momentum transfer 0 

4 Conclusions [TUI 


1. Introduction 

Stimulated scattering of light or particles from condensed-matter systems - solids, 
liquids, and gases - is a powerful tool for providing fundamental insight into the structure 
of matter. Elastic scattering of x-ray photons has permitted to disclose the atomic order 
and electron distribution in crystalline solids, as well as the arrangement of atoms in 
molecules [T]. Similarly, inelastic neutron scattering has unveiled the phonon spectrum 
of superconductors and the superfluidity of liquid helium [2]. 

In cold atomic systems, inelastic scattering of photons - also known as Bragg 
spectroscopy - has been used to study harmonically trapped three-dimensional (3D) 
Bose-Einstein condensates (BECs) [21 IH El [6], strongly interacting BECs across 
a Feshbach resonance [7], and strongly interacting fermions El i, through direct 
observation of the net momentum imparted to the system. The transferred momentum 
is easily measured in this kind of settings, since the atomic density distribution, observed 
after time-of-flight in the far-held regime, directly rehects the in-trap momentum 
distribution. Instead, cold atoms in optical lattices have been investigated by measuring 
the increase of energy following a modulation of the lattice amplitude where the 
excitation has zero momentum [10], and with scattering experiments where the 
excitation has non-zero momentum [mini OS). The energy of a condensate, even 
in the presence of shallow lattices, is easily extracted from the time-of-hight density 
distribution of the gas [TT], whereas the energy of strongly-interacting systems realized 
in deep optical lattices - as a Mott insulator - is not directly accessible, unless with 
single-site resolution experiments na. In the case of deep optical lattices, the energy 
excess produced by the Bragg perturbation can be measured by lowering the lattice 
depth, i.e., driving the system in a less interacting regime, and let it thermalize mm- 

In the linear response regime, both energy and momentum transfer are related to 
the dynamic structure factor [Tn[, which carries key information about the dynamical 
behaviour and correlations of the system. However in trapped condensates, while the 
energy is a conserved quantity, momentum is not conserved due to the presence of the 
trap. Thus, if the Bragg pulse duration is not negligible compared to the inverse of the 
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trap frequency, the momentum imparted by the Bragg beams can be affected by the in¬ 
trap dynamics USE], complicating the connection between the momentum transferred 
and the dynamical structure factor. On the other hand, a short Bragg pulse would result 
in a limited spectral resolution. 

In this work, we use inelastic light scattering for accessing the dynamical structure 
factor of an array of one-dimensional (ID) Bose gases. The dimensionality of the system 
plays a crucial role: in ID quantum systems, correlations - which directly reflect on 
the dynamical structure factor - bring to peculiar phenomena, such as fermionization 
of strongly-interacting bosons [191 ISHl El], or spin-charge separation of interacting 
fermions [22], which do not have any higher-dimensional equivalent. Moreover, for ID 
systems, the mechanisms and characteristic times of thermalization are currently under 
investigation [221 El] • This may affect the measurement of the energy transferred via 
light scattering. On the other side, momentum measurement of ID trapped gases may 
be influenced by the in-trap dynamics, as above mentioned. The purpose of this work 
is to investigate in experiment the relation between energy and momentum imparted to 
an array of ID gases due to Bragg scattering, in a typical regime of parameters mm, 
and to discuss the effect of the in-trap dynamics on the transferred momentum. 

The paper is organized as follows. In Sec. we focus on the comparison of the 
response of the array of ID gases to the scattering experiments in terms of energy 
deposited and momentum boost imparted to the system. We present the experimental 
setup and discuss the results, obtained in a regime of weak perturbation that is well 
described in the framework of the linear response theory. We also directly compare the 
susceptibility of this system to that one of a 3D non-interacting condensate. In Sec. 
we study the effect of the in-trap dynamics on the momentum transfer by recording the 
evolution of the response of the system in time, after the Bragg excitation. 

2. Energy and momentum transferred to an array of ID gases 

2.1. Experimental setup 

We produce an array of ID gases by loading a Bose-Einstein condensate of ®^Rb in a 
two-dimensional optical lattice created by two mutually orthogonal laser standing waves 
of wavelengths Xl = 765 nm. The loading is performed with an exponential ramp of 
tr =250 ms, with time constant fj,/3. The hnal depth of the lattice is Vl = 30Er, with 
Ey. = h‘^/{2mX\), m being the atomic mass and Xl the lattice wavelength. This value 
is chosen to be high enough to freeze the transverse degree of freedom of each ID gases 
(the radial trapping frequency is u}± = 2 ti x (42 ± 2) kHz), and suppress the tunneling 
of particles between different tubes on the timescale of the experiment. 

The equilibrium state of the system is completely described by two dimensionless 
parameters [2S|, that is, (i) the interaction parameter 7 = rngiD/{h?p), where giD 
is the one-dimensional interaction strength [ 26 ], and (ii) the reduced temperature 
t = 2mkBT/{hpY, which depend on density p and temperature T. 
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The significant parameters that characterize the array can be estimated by rescaling 
the interparticle interaction strength as in Ref. [13]. We estimate the array to consist 
of about 4 X 10^ ID tubes, and the central tube to have 7 ~ 1 and f ~ 1, with density 
p ~ 5 pm“^ and chemical potential /i/h ~ 3.6 kHz. 

The study of this system is carried out by imparting a perturbation to the array 
of ID tubes given by two simultaneous off-resonance laser pulses with time duration 
= 3 ms, which determines an interaction-time broadening of ~ 150 Hz. Note that 
tB — T/4, with T = 27i/oJx being the trap period along the axis of the tubes. The laser 
light is detuned by 200 GHz from the ®^Rb resonance. The two beams have tunable 
relative detuning Aco (up to tens of kHz), and produce a moving Bragg grating with 
amplitude Vb = h x 900 Hz. The wavevector of the Bragg grating is adjusted to be 
along the axial direction of the tubes, and it is hxed at g = (7.3 ± 0.2) pm“^. In the 
experiments, we vary the energy hu of the excitation by tuning Au (being u = |Aa;|) 
and we measure the response of the system, in terms of energy and momentum transfer. 


2.2. Results 

For measuring the energy transfer, after the Bragg pulse we lower the lattice height to 
Vl = 5Er, where the whole gas is expected to thermalize in a superfluid regime. After 



|A(o|/(27r) (kHz) 


|A(o|/(27i) (kHz) 


Figure 1. Energy transfer, (a) Experimental timing used for measuring the energy 
spectrum. The optical lattice is raised to 30Er in 250 ms. During the holding time 
the Bragg beams are shined onto the atoms, then the whole system is let thermalize 
by lowering the lattice depth to 5Er. After 5ms the trap is released, and the energy 
is measured from the time-of-flight density distribution, (b) Energy transferred via 
Bragg scattering to an array of ID gases in an optical lattice of height Vl = 30 Er, as 
a function of the excitation frequency. The signal is obtained from the squared width 
of the central peak of the time-of-flight density distribution -|- ^ay). The fitting 
curve is given by a gaussian function multiplied by the frequency wGioj). (c) Increase 
of the rms size along the axis of the tubes {aR) and along a perpendicular axis (cy), 
as a function of the excitation frequency. 
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5ms the gas is released and let expand ballistically for a time-of-flight ttof = 25ms, 
then we record the density distribntion of the atomic cloud. The experimental timing is 
sketched in the Fig. [^a). From the time-of-flight images, we extract the squared width 
of the central peak of the resulting interferogram cr^ = cr^ -|- 2(7^ [27] and subtract the 
background cTq, corresponding to the value measured in the absence of the Bragg pulse, 
in order to obtain the experimental signal. This quantity is proportional to the energy 
imparted to the system, as previously demonstrated [13] . The latter, in turns, is related 
to the dynamical structure factor S{q,u) of the system through the relation [2H] 

AE{q,u)= tBU}S{q,u), (1) 

valid in the linear response regime. The measured energy spectrum, normalized to its 
integral, is shown in Fig. [^b). In order to verify that the system has thermalized after 
the Bragg pulse, in Fig. [^c) we plot separately the increase of the rms size observed in 
each direction {Aax and Aay): in spite of the symmetry breaking induced by the Bragg 
perturbation imparted along the axis of the ID gases (x direction), Aa^ and Acr^^ show 
the same dependence on the Bragg frequency, indicating that an efficient thermalization 
process has occurred during the ramping down of the lattices. 

In the experiment, we also measure the total momentum imparted to the same 
system. To this purpose, after the Bragg pulse, the atoms are abruptly released directly 
from the trap, as represented in Fig. [^a), so that in-trap momentum is mapped into 
the atomic density distribution after time-of-flight. When the Bragg perturbation is 
on resonance, momentum is transferred efficiently, and the time-of-flight images of the 
density distribution exhibit a small cloud of excited atoms ejected from the main cloud. 
Figure l^b) shows the evolution of the normalized density profiles n{x) integrated along 
the line of sight and along the y direction (orthogonal to the axis along which momentum 
is transferred) with the relative detuning between the Bragg beams. From the time-of- 
flight images, the net moment boost AP{u, q) is obtained by measuring the displacement 
of the center of mass - relative to the unperturbed position - as [29] 

TTl f 

AP{q, u) = -— X (n(x, y) - no{x, y)) dxdy , (2) 

ttof J 

where n{x,y) and no{x,y) are the density profiles integrated along the line of sight 
and normalized to the unity, with and without the Bragg excitation, respectively. The 
experimental spectra normalized to the momentum of the excitation hq, AP{q,uj)/{hq), 
are shown in Fig. [^c). Filled (empty) dots correspond to positive (negative) 
momentum, along the axis of the tubes. 

As remarked in Refs. [ig [m ED], momentum is a conserved quantity only in the 
absence of any external trapping potential. In the latter case, for a perturbation in 
the linear response regime and with ojIb 3> 1, AP{q,uj) is related to the dynamical 
structure factor through the following relation 

AP{q,u)= tBqS{q,u). (3) 
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Figure 2. Momentum transfer, (a) Experimental timing. The optical lattices are 
abruptly switched off immediately after the Bragg pulse, and the momentum transfer 
is measured from the center-of-mass shift of the whole cloud, (b) Normalized density 
profile n{x) along the axis of the ID gases - integrated along the line of sight z and the 
y direction - for different values of the relative detuning of the Bragg beams Aw/(27r) 
(false colours). On resonance, the central peak is depleted and the atoms are ejected 
along the x direction, (c) Momentum transferred to an array of ID gases. The filled 
(empty) dots, at positive (negative) momentum along the axis of the tubes, correspond 
to Aw < 0 (Aw > 0). The Bragg parameters of this measurement are the same as 
described for the measurement of the energy transfer. 


For a trapped gas with axial trapping frequency Ux significantly smaller than the radial 
one, this equation still holds in a wide range of parameters, provided that S> 1 and 
UxtB < 1 In the present case, the hrst condition is well satished as utB — 80 on 
resonance, while we have u:xtB — which does not satisfy the second condition. Thus 
the comparison between the quantities extracted from the measurements of energy and 
momentum transfer is not straightforward. In order to address quantitatively this issue, 
we £t AP{q,u) with a gaussian function G(c(;), where center, width and amplitude are 
free parameters, and AE(q, u) with oo G{u), as follows from Eq. The gaussian centers 
obtained from the measured spectra are respectively (4.5 ± 0.2) kHz for the momentum 
transfer, and (4.3 ± 0.3) kHz for the energy transfer, with their widths being (2.5 ± 0.2) 
kHz and (2.3 ±0.2) kHz, respectively. These results are consistent within the error bars, 
allowing us to conclude that, with this choice of parameters, both these experimental 
approaches measure the same quantity. 

2.3. Comparison with the response of a 3D condensate 

As a reference system, we also measured the transferred momentum of a 3D expanded 
condensate, since its response is well described by a non-interacting model. The 
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interparticle interactions are suppressed indeed by letting the BEC free fall for 5 ms 
of time-of-flight before performing the scattering experiment. For direct comparison, in 
Fig. [^we show the experimental spectrum of the array of ID gases obtained as previously 
described (Fig. [^a)), and the spectrum of the 3D non interacting condensate (Fig. 
[^b)). In both Figs, [^a) and 3(b), the signal is normalized to the Bragg strength Vits¬ 
in this Figure, we also report the exact solution for a free-particle system (continuous red 
line), which shows an excellent agreement with the experimental data. This prediction is 
obtained by solving the Schrodinger equation BD in the presence of the Bragg potential 
V{x, t) = 6{t—tB)VB cos{qx—ujt) (for f > 0), and does not contain any fitting parameter. 
From the comparison of the two spectra, we can notice that the response of the ID 
tubes is much broader, as a consequence of interactions [13]. Moreover, the comparison 
between the two amplitudes indicates that the susceptibility of the ID tubes is about 
35 times lower than the one of the 3D non-interacting condensate. 

The very low susceptibility of the array of ID tubes, with the Bragg parameters 
that we have used, is a first indication that the response of the system lies in the linear 
regime of weak perturbation, where the relations in Eq. ([^ and Eq. (|^ are expected 
to hold. We have also verified the behaviour of the experimental signal as a function of 
Bragg power in a range that includes the used value, and demonstrate the dependence 
of the signal on Vb to be quadratic, as expected in the framework of the linear response 
theory [T6] . 


3. Effect of in-trap dynamics on the momentum transfer 


As previously discussed, momentum is not a conserved quantity in the presence of a 
trapping potential. Therefore, the measurement of the momentum following a Bragg 
excitation can be affected by the in-trap dynamics before the release [I7l|30|. For the 
cases discussed so far, the momentum transferred has been measured immediately after 
the Bragg pulse, see Fig. For the ID gases, even if ~ T/4, as seen in Sec, 2.2, AP 
and AP carry the same information, then we can infer that no appreciable effect of the 
dynamics during the Bragg pulse has been observed. 

Now, we characterize more in depth the effect of the in-trap dynamics on the 
response of the ID gases. To this purpose, we measure AP{q,u)/{hq) at variable time 
after the end of the Bragg perturbation. In Fig. |^a) we can observe a modification of 
the system response with time. The Bragg time-duration is fixed to the value = 3 
ms, and the total holding-time of the atoms in the lattice trap {tn) is kept constant, 
while we vary the time tc between the end of the Bragg pulse and the release of the 
trap from 0 ms up to 10 ms, as sketched in the inset of the figure. During the first 2 ms 
after the Bragg pulse, the total spectral weight of the signal undergoes a suppression 
with time (note that the vertical scale is the same in all the panels of the Figure), which 
eventually results in a negative amplitude at tc = 10 ms. Remarkably, the shape of the 
signal at fc = 4 ms is asymmetric and qualitatively different from the other cases. 

The latter behavior can be understood by considering the effect of the in-trap 
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Figure 3. Comparison between the transferred momentum AP to (a) an 
array of trapped ID gases and (b) a non-interacting 3D condensate. The 

horizontal scale represents the frequency difference between the Bragg beams. The 
spectrum of the array of ID gases has been obtained using a Bragg parameters ts = 3 
ms and Vb = hx 900 Hz, whereas the spectrum of the 3D condensate has been obtained 
with tB = 0.5 ms, and Vb = hx 540 Hz. The amplitude of both the spectra has been 
rescaled by the pulse strength Hjfs to directly compare the susceptibility of the two 
systems. Note that the vertical scale in the first graph is 10 times smaller than in 
the second one. A fit of the experimental spectrum of the ID gases with a sum of 
two gaussian functions (blue curve) is shown in Fig. (a) and reported also in Fig. 

(b) for highlighting the comparison between the response of the two systems. The 
red continuous line in Fig. (b) is the solution of the time-dependent Schrodinger 
equation for a non-interacting gas, given the value of Vb/ h, with no free parameters 
(see text). 


dynamics during the Bragg pulse. Let us consider a system of N interacting particles 
trapped in a harmonic potential, described by the following Hamiltonian 

N 

i=l 

The evolution of the total momentum and position operators along each spatial 
directions can be easily obtained from the Heisenberg equations as = (—i/h)[iL,pa] = 
mUj'^Xa and Xa = {—i/h)[H,Xa] = —pa/m, with Oa = = 1)2,3, O = x,p). 

For the hrst relation we have used the fact that = —dxi^VpN. Then, restricting 

the discussion to the ID case, it is straightforward to get that the average momentum 
evolves in the trap as 

{p){t) = -mUx{x)osm{uxt) + m{x)ocos{uxt) (5) 

where (T)o and (x)o are the average values of the density and velocity distributions at 
time t = 0"^ immediately after the end of the Bragg pulse. We remark that this result is 
valid in general for any interacting system, regardless of the temperature, the statistics 
(being the particles bosons or fermions) and the dimensionality of the system. In fact, 
it is a well known result that the dynamics of the center of mass in the presence of 


j<i 


X 


-Xj) 


(4) 
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Figure 4. In-trap dynamics after the Bragg excitation. AP/{hq) is reported 
as a function of |Aa;|/(27r), for different holding times in the trap after the Bragg 
excitation: tc = 0,1,2,4,10 ms. The continuous blue lines are guides for the eyes. 
The curve corresponding to t = 0 ms is also reported in the other panels (as a dashed- 
dotted line) for comparing the signal shapes and amplitudes. Inset: experimental time 
sequence. The total holding time in the lattices is kept constant {tn = 30 ms). 


harmonic trapping is decoupled from the internal degrees of freedom of the system (see 

e.g. [H]). 

Let us now turn to the effect of the Bragg pulse, that we assume of the form 
g(t)VB cos{qx — ojt). First, let us consider an instantaneous Bragg pulse described by 
g{t) = 6{t). After the pulse, at f = 0+, the density distribution is basically unperturbed 
((x}o = 0). Then, the Bragg perturbation only affects the initial velocity distribution 
n(xo), so that its mean value is (x}o = {Nb{oj) /N)hq/m, where Nb/N is the ratio 
of the number of diffracted atoms to the total number of atoms, and depends on the 
excitation frequency oj/{2'k). As follows from Eq. (|^, in this case {p){t) vanishes exactly 
for f = T/4, T being the period of the trap, for any excitation frequency. 
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Instead, for a finite dnration of the Bragg pnlse (and in particnlar, if Ib is 
comparable with the trap period), even the spatial distribntion of the atomic ensemble 
may nndergo modihcations dnring the Bragg pertnrbation, depending on the excitation 
freqnency. This makes the initial valne of the center-of-mass {x)o in Eq. (|^ non 
vanishing and ca-dependent, therefore affecting the following dynamics and changing 
the shape of the signal. 

As an example, let ns consider the simple case of a Bose-Einstein condensate in a 
single, qnasi one-dimensional tnbe, in the mean-held regime. In this case, the response 
of the system to the Bragg pnlse can be easily obtained by solving the following ID 
Gross-Pitaevskii eqnation (t > 0) 


ihdt'ip 



+ Vho{,x) + 6{t - tB)VB cos{qx 




V’, (6) 


where qid = g/{27ra±), g = 4nh‘^a/m being the 3D interaction strength, a the scattering 
length for ®'^Rb, and a± = \/h/{muj±) the oscillator length in the transverse directions. 
In this specihc example we consider Vb/Ii = 120 Hz, = 3 ms, = 2ti x 60 Hz, 
a;_L = 27 r X 42 kHz, and an array of tnbes that corresponds to the typical experimental 
conhgnration. The response of the system at different evolntion times in the trap 
is shown in Fig. where the meanheld predictions are also compared to the non¬ 
interacting case [32]. This hgnre shows that indeed, as follows from Eq. (|^, the 
response patter at t = T/2 is reversed with respect to that at t = 0, the evolntion 
being periodic in time. For intermediate times, the shape of the signal is non trivial, 
depending on the relative weight and on the specihc shape of the transferred momentnm 
and the center of mass position as a fnnction of the Bragg freqnency |Aa;|, at t = 0. In 
the non interacting case, (p)o(|Aci;|) and (x)o(|Aci;|) are centered at the same valne and 
almost symmetric aronnd that point, so that the same symmetry property is preserved 
dnring the evolntion. Instead, the response of an interacting condensate is characterized 
by a distribntion of the center-of-mass position that is peaked at higher freqnency 
with respect to the corresponding transferred momentnm, and this affects dramatically 
the shape at intermediate times. In particnlar, at t ~ 0.15T the signal shape has a 
characteristic sinnsoidal-Iike form, whereas at f = 0.25T it corresponds exactly to the 
reverse of the initial center-of-mass distribntion {x)o(|Aa;|). Note that the signal at 
t ~ 0.15T is analogons to that observed experimentally at t = 4 ms, see Fig. thongh 
in that case the signal hrst reverses in the low freqnency region. We attribnte this 
difference to the fact that the experiment is performed in the strong-conpling regime, so 
that the response of the system to the Bragg pertnrbation is expected to be snbstantially 
different (we remark that a precise simnlation of the dynamics of strongly correlated ID 
systems nnder the effect of a Bragg pertnrbation can be very demanding, see e.g. [33]). 


4. Conclusions 

In conclnsion, we have investigated the response of an array of ID gases, comparing 
energy and momentnm transfer in Bragg spectroscopy experiments. In the presence of 
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Figure 5. In-trap dynamics of a condensate in the mean field regime. 

Theoretical prediction for the value of AP/{tiq) (blue squares) as a function of the 
Bragg frequency for an array of quasi ID BECs, for different holding times in the 
trap after the Bragg excitation: t = 0, 0.15, 0.25, 0.5r (T being the trap period). The 
dotted-dashed (red) lines in the top panel represent the initial distributions of the 
center-of-mass (x)o(|Aa;|). The same line is shown inverted aX t = T/4 (see text). The 
left and right columns correspond to (a) the non interacting and (b) interacting case 
in the mean-field regime, respectively. 


an external trapping potential along the axis of the tnbes, even if increasing the pnlse 
time enhances the spectral resolntion, the presence of the trap in principle provides an 
npper limit to the pnlse dnration. In addition, onr experiment reveals that, in a regime 
of parameters well described by the linear response theory, and for time-dnration of the 
Bragg pertnrbation smaller than a qnarter of the trapping period, the proportionality 
relation between the momentnm transfer and the dynamical structnre factor is well 
respected. Moreover, we show that the in-trap dynamics during the Bragg pulse affects 
noticeably the response of the system. This analysis can be useful for interpreting the 
results of scattering experiments also in other more complex settings of ultracold gases 
in optical lattices or disordered potentials. 
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